Berry Phases of Classical Trajectories in the Presence of Hard- Wall Boundaries 
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We study the quantum propagator in the semiclassical limit with hard-wall potentials. We show 
that, upon each reflection by the hard wall, a Berry phase vr is accumulated and leads to interfer- 
ences between different classical trajectories. Including the Berry phase caused by the hard walls, 
the modified Van Vleck's formula is derived. We also discuss the close relations among quantum 
statistics, discrete gauge symmetry, and hard-wall constraints. Most of all, we formulate a new quan- 
tization rule that applies to both smooth and sharp boundary potentials. It provides an easy way 
to compute quantized energies in the semiclassical limit and is extremely useful for many physical 
systems. 



I. INTRODUCTION 

Path integral provides an alternative approach to for- 
mulate quantum mechanics The quantum prop- 
agator G{x,x';T), that is the key quantity in quantum 
mechanics, is shown to equal the summation over all pos- 
sible paths with the same end points. In the semiclassical 
limit, the dominant contribution comes from classical tra- 
jectories and fluctuations around them [^-^. Within the 
stationary phase approximation including fluctuations up 
to the quadratic order, the quantum propagator can be 
approximated by the Van Vleck's formula . In general, 
there would be many classical trajectories that satisfies 
the same boundary conditions, and the Berry phase in- 
terferences between them are important . 

The interference effects among the classical trajectories 
were not included until Gutzwiller's work that enables 
the earlier short-time result of Van Vleck to be extended 
past the first conjugate point [||. By Morse's theorem, 
the second variation, considered as quadratic fluctuations 
around a given trajectory from x' to x in time T, has as 
many negative eigenvalues as there are conjugate (turn- 
ing) points along the trajectory. These conjugate points 
give rise to a Berry phase j/7r/2 for the trajectory, where 
1/ is the total number of conjugate points along the tra- 
jectory, or sometimes referred to as the Maslov or Morse 
index [0. 

Not only elucidating the crossover between classical 
and quantum mechanics, the semiclassical limit also pro- 
vides a convenient way to calculate the bound state en- 
ergy. Instead of solving the Schrodinger equation di- 
rectly, the bound state energies can also be computed 
by the WKB approximation |jll[]. However, the method 
is only applicable when the confinement potential is rea- 
sonably "smooth" (compared with the relevant energy 
gradient we are interested in). In this paper, we study 
the semiclassical limit in the presence of hard-wall po- 
tentials that is beyond the validity of the WKB approxi- 
mation. We found that upon each reflection by the hard 
wall, a Berry phase tt accumulates and eventually leads 
to an additional phase rn, where r is the number of re- 
flection points. In order to account for the interference 



effects among classical trajectories correctly, we rederive 
Van Vleck's formula with an extra Berry phase correc- 
tion due to hard- wall boundaries. An elegant proof by 
mirror projection explains the Berry phase tt and reveals 
the close connections among quantum statistics, discrete 
gauge symmetry, and the hard-wall boundary. 

Following the standard stationary phase approxima- 
tion and making a Legendre transformation of the time 
variable in the quantum propagator to the energy vari- 
able, we are able to generalize the Einstein-Brillouin- 
Keller (EBK) quantization rule |T^-[l^ with a hard-wall 
correction term 



j) ^/2m[E - V{x)]dx = 27r(n - 



4 + 2^' 



(1.1) 



where v is the number of conjugate points and r is the 
number of reflection points. The usual WKB approxima- 
tion is the special case = 2 a nd r = 0. The modified 
EBK quantization rule in Eq. (1.1) relaxes the require- 
ment of the potential smoothness in the WKB approxi- 
mation. If the trajectory bounces from a smooth poten- 
tial, it is counted as a conjugate (turning) point. If the 
trajectory gets reflected by a hard-wall like boundary, it 
is counted as a reflection point. This is of great advantage 
because many physical systems including quantum dots, 
quantum wells. Hall bars, electronic wave guides, etc., 
have both hard-wall-like potentials (from sample edges) 
as well as smooth potentials (by applying external fields) 
at the same time. 

The paper is organized in the following way. In Sec. II, 
we introduce the Van Vleck's formula and apply it to two 
simple systems without and with hard-wall boundaries. 
We explicitly show that the Van Vleck's approximation 
is incorrect in the presence of hard walls. In Sec III, 
We show that the Berry phase due to the refiection by 
a single hard-wall boundary is tt and derive the modified 
Van Vleck's formula. In Sec. IV, a more elegant proof 
by mirror projection is presented. Finally, in Sec V, we 
derive the key result of this paper - the modified EBK 
quantization rule. We apply it to physical systems with 
both smooth and hard confinement potentials and show 
that the modified term is necessary to obtain the correct 
energy levels. Then a brief conclusion follows. 
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II. QUANTUM PROPAGATOR AND CLASSICAL 
TRAJECTORIES 



G(x, x; T) 



2mT ■ 



^exp [z^(a;-x' + nL)2j (2.5) 



In the path integral formahsm [Q, the quantum prop- 
agator equals the sum over all possible paths with the 
same end points, 



G{x,x'\T) 



'D[x\ exp 



I I L(x,x,t)dt 



(2.1) 



where the measure T^lx] denotes all possible paths with 
end points x{0) — x' and x{T) — x. In the semiclas- 
sical limit, the phase inside the path integral oscillates 
rapidly except in the neighborhood of the classical tra- 
jectories. Within the stationary phase approximation in- 
cluding fluctuations up to quadratic order, the propaga- 
tor is approximated by the Van Vleck's formula. 



G(x,x';T) 



^^/C^exp[iAp- ivp-], (2.2) 



where Ap{x, x'; T) is the action of the classical trajectory 
starting from a;(0) — x' and ending at x{T) = x, and the 
subscript p denotes all classical paths with the desired 
end points. The strength of the quadratic fluctuations 
PI around the classical trajectory is 



Cri 



dxdx' 



(2.3) 



Finally, the total number of conjugate (or turning) points 
along the classical trajectory is denoted by v. Notice 
that, for each conjugate point, there is a 7r/2 Berry phase 
associated with it. Van Vleck's formula provides a com- 
pletely classical approximation of the quantum propaga- 
tor, in the sense that all relevant elements can be com- 
puted from the classical trajectories. 

A straightforward example of the Van Vleck's formula 
is a free particle moving on the a finite ring with length 
L. There are infinite classical paths which satisfy the 
conditions x(0) = x' and x{T) — x. The total (route) 
distance of each classical trajectory ' -t- riL, 

where n is an integer. The action for each trajectory is 



An{x,x ;T) — i^i^ ' 



nLf 



(2.4) 



Taking the derivative of the action, the strength of fluctu- 
ations around each trajectory C„ = m/T is independent 
of the end points and the choice of trajectories. Since 
the particle moves at constant velocity, it is obvious that 
there is no conjugate point along any classical trajec- 
tory and thus Un = 0. Besides, because the fluctuations 
of the classical trajectory of a free particle are exactly 
quadratic, we expect the Van Vleck's formula to be ex- 
act for this system. 



This infinite sum can be re- written in terms of its Fourier 
function with the use of Poisson summation formula in 
Appendix A. Notice that 



.f{y) 



^ F{p) = . &-»feV4a+^fc;3 (2.6) 

V a 



Choosing a = L, the summation over coordinate y = na 
can be turned into the summation over momentum A;„ = 
2m: /L. The propagator is then 

G{x, x;T)^^Yl ^Miknix - x') - iE^T], (2.7) 



where fc„ = 2m: /L is the quantized momentum and 
En = fc^/2m is the quantized energy. It is obvious that 
the propagator G{x,x';T) calculated by the Van Vleck's 
formula is exact in this case. 

Let us now apply the Van Vleck's formula to another 
physical system ~ a free particle bouncing back and forth 
between two hard walls. We would calculate the propa- 
gator explicitly and show that the Van Vleck's formula 
leads to incorrect results. 



FIG. 1. Classical trajectories in the presence of two hard 
walls. On the left is a trajectory with even reflection points 
r — 2, while the right with odd reflection points r — 3. 

The trajectories in this problem can be classified by 
the number of collisions with the hard walls, as seen in 
Fig. 1. For those trajectories that collide with the hard 
walls even times, the route distance is = x — x' + 2nL, 
while the distance is — x + x' + 2nL for trajectories 
that collide with the walls odd times. The action for each 
trajectory can be computed straightforwardly 

Al{x,x'-T) = —{x-x' + 2nLf, (2.8) 
Al{x,x'-T)^^{x + x' + 2nLf. (2.9) 

Here A'^/°{x,x' ;T) denotes the action for trajectories 
with even/odd reflection points. The fluctuations along 
all trajectories contribute the same C„ = m/T as in the 
previous example. For an one-dimensional motion, a con- 
jugate point is identifled as the position where the veloc- 
ity vanishes. However, for a free particle bouncing back 
and forth between two hard walls, the velocity is constant 
up to a minus sign and does not vanish at any point along 
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the classical trajectory. Thus, the number of conjugate 
points is zero, v = 0. 

The propagator without any Berry phase interference 

is 



Gvvix, x'- T) = { ^Mi^{x -x' + 2nLf] 



+ exp[i — (a; + a;' + 2ni)2]|. (2.10) 

Both infinite sums can be turned into summations over 
discrete momentum again by mean of the Poisson sum- 
mation formula. The prefactors cancel as in the previous 
example and we are left with the simple result, 

^ oo 

Gvv{x,x'\T) = - ^exp[-z£;„t]x 

n=0 

X I cos[fc„(a; — a;')] + cos[A;„(a; + a;')]|, (2.11) 

where fc„ = mi/ L is the quantized momentum and £"„ = 
kf^/2m is the quantized energy. Combining two cosines 
would leads to cos(A:„a;) cos(fc„a;'), while the correct form 
should be sin(/c„x) sin(fc„a;'). Indeed one can recover the 
exact answer (with all prefacto rs rig ht!) if we change the 
sign of the second term in Eq. ( 2.11 ). That is, only if we 
assign an extra Berry phase tt to trajectories with odd 
reflection points, will the modified Van Vleck's formula 
become correct! 

In the following section, we study the path integral 
formalism in the presence of a single hard-wall boundary 
and show that there is an additional Berry phase. 



III. BERRY PHASES DUE TO HARD WALLS 

Consider a particle moving under the influence of a 
regular potential V(x) and a hard- wall potential Vc{x). 
The Hamiltonian is 



H = ^ + V{x) + V,{x), 
where Vc{x) is the hard- wall potential at a; = 0, 



Vcix) 



0, a;>0; 
CO, X < 0. 



(3.1) 



(3.2) 



The regular potential is treated in the ordinary way while 
the hard-wall one is viewed as the depletion of Hilbert 
space. The complete set of the Hilbert space is now re- 
duced, 



dr\r) (r| — 1, 



(3.3) 
(3.4) 



It would become clear later that the phase 4> is associ- 
ated with the Berry phase in the path integral. Slicing 
the time interval T into infinitesimal pieces and inserting 
complete sets of the coordinate space, the propagator is 



G{r,r';T) = {r\e-^"^y: 



.oo N-l 

I dr^ TT (r„+i|e 



(3.5) 



where rj^ = r and tq = r' are all positive. Each ma- 
trix element in the product is computed by inserting the 
complete set in momentum space into Eq.(3.4), 

dpn 

Jp,t{r„+i-x„)-itj>„ 



{rn+i\e-"-"\r„ 



= / exp [-zei/„] 



E 



X > e^' 



(3.6) 



where the phase = for Xn = Vn, and (j) = tt when 
Xn = — r„. Since Xn = i?'™, the two terms can be com- 
bined and lead to the unconstraint integral over a;„. After 
changing the constrained variable r„ to Xn, it is conve- 
nient to write the Berry phase 0„ in the following way 



(/)„ = 7r[e(x„+i) - 6(x„)]. 



(3.7) 



Notice that the Berry phase is zero if the path does not 
pass through a; = in the infinitesimal time interval dtn 
and TT if the path passes through. The integral over mo- 
mentum can be carried out easily and the propagator is 



G{r,r';T)= E e*'^'' / exp[iyl(r, x'; T)]. 

x'=±r' 



(3.8) 



The total Berry phase (j)B = Tr[Q{r)—Q{x')] is a boundary 
term and can be pulled out of the path integral |^ . The 
paths are divided into two topologically distinct classes. 
For all possible paths starting from r to r', the Berry 
phase is zero, while for those starting from r to — r', the 
Berry phase is tt that causes a minus sign. The classical 
trajectories among the paths can be then classified in the 
same way. Furthermore, trajectories with end points r 
and r' can be identified as trajectories (in the physical 
half plane) with even reflection points and those with 
end points r and — r' are trajectories with odd reflection 
points. 

Therefore, in the semiclassical limit, the Van Vleck's 
formula is modifled with an extra Berry phase term, 

G(r, r'; T) ~ E V^expMp 

yzTTi — 
p 



0=0, TT 



(3.9) 

where Tp is the number of reflection points. The proof for 
more than one hard wall is straightforward but tedious. 
A more general and elegant proof would be given by mir- 
ror projection in the next section. However, I would like 
to emphasize that the proof given in this section explains 
clearly the origin of the minus sign ~ an extra Berry phase 
due to hard-wall boundaries. 
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IV. MIRROR PROJECTION 

In the previous section, we treat the hard-waU bound- 
ary as depletion of the Hilbert space. An alternative way 
is to view it as a discrete Zi gauge symmetry of the wave 
function 

i/,(a;) = -i}{-x). (4.1) 

The minus sign is chosen here to make the wave func- 
tion vanishes at x = so that the boundary condition 
■0(0) = is always satisfied. Since the propagator can 
be written down as the summation of eigenfunctions, 
G{x,x';T) = V'«(a;)'(/'^(2;')exp[-i£'„T], where V«(a;) 
is the eigenfunction with eigenenergy i?„. The discrete 
gauge symmetry of the wave function implies that the 
quantum propagator has the symmetry 

Gix, x'; T) = -G{x, -x'; T). (4.2) 

Now choose both x = r and x' = r' to be positive, 
the propagator can also be viewed as the wave func- 
tion G{r,r';T) = ij}r'{r,t) that satisfies the Schrodinger 
equation with a delta function source at {x,t) = {r',0). 
The propagator Go{r, r'; T) without the hard- wall bound- 
ary satisfies exactly the same differential equation except 
that the boundary condition at a: = is not met. Notice 
that the mirrored propagator Go(r, r'; T) = Go{r, — r'; T) 
satisfies the Schrodinger equation without the source 
term since the delta function S{r + r') ~ for positive 
coordinates. Therefore, the propagator that satisfies the 
correct boundary condition is constructed as 

G(x, x'; T) = Go(x, x'; T) ~ Go{x, x'; T). (4.3) 

The above result is equivalent to Eq. ( b.Sh. I t is obvious 
that the discrete gauge symmetry in Eq. ( [4.2[ ) is satisfied. 
This method is just the familiar mirror charge trick in the 
classical electromagnetism. The generalization now be- 
comes clear. In the presence of more than one hard walls, 
the propagator is just the sum of all path integrals from 
all mirror points of x' to the final x, and the Berry phase 
correction is either or tt depending on how many mir- 
ror projections are needed to reach the particular mirror 
point of x' . In the example of two hard walls, there are in- 
finite mirror points. Each mirror projection corresponds 
to a reflection of the classical trajectory. Thus, for those 
trajectories that are reflected even times, the Berry phase 
is simply zero, while for those with odd reflection points, 
the Berry phase is tt and results in a relative minus sign. 

Since we can solve the hard-wall boundary by dis- 
crete gauge symmetry, we might as well go the other 
way around. It is possible to replace the quantum statis- 
tics between particles by the hard-wall boundaries. Let 
us consider the simplest case - two interacting particles 
with either bosonic or fermionic statistics. The discrete 
gauge redundancy is 

V'(x) = e*'^V(-a;), (4.4) 



where = for bosons and tt for fermions. The dis- 
crete gauge symmetry is removed by imposing a hard wall 
xi — X2 in the configuration space, and a Berry phase </> 
accumulates upon each reflection due to the hard wall. 




(a) (b) 
FIG. 2. Classical trajectories of two particles whose quan- 
tum statistics is replaced by the equivalent hard wall at 
xi = X2 in the configuration space. In part (a), a direct 
trajectory is shown and, in part (b), the shown reflected tra- 
jectory is equivalent to exchanging two particles which results 
in an extra Berry phase. 

Classical trajectories are classified into two categories 
- the direct path and the reflected one as shown in Figure 
2. If we extend the reflected trajectory into the unphysi- 
cal regime inside the hard wall, as shown in Fig. 2(b), the 
reflected trajectory is equivalent to an exchange between 
two particles. This approach would be useful when study- 
ing few interacting quantum particles, e.g., two strongly 
interacting bosons or fermions bouncing back and forth 
between two hard walls. In the semiclassical limit, we 
can safely ignore the quantum statistics by solving all 
classical trajectories inside a specific triangle in the two- 
dimensional configuration space. 

V. MODIFIED EBK QUANTIZATION RULE 

The most powerful use of Van Vleck's formula is that 
it leads to the EBK quantization rule in the semiclassical 
limit. One notices that, if we set x = x' in the propagator 
and integrate over all possible x, it results in the quan- 
tum partition function Z(T) = '^^e:xjp[—iEnT]. The 
energy levels can then be identified as the singularities of 
Z{uj) which is the Fourier transformation of the partition 
function. Within stationary phase approximation, it can 
be shown that the total Berry phase § pdq — ivi: J2 (in 
the absence of hard- wall boundaries) is quantized ||] and 
leads to the EBK quantization rule, 

/p<„ = 2.(„+^,, (5.1) 

where v is the number of conjugate points along the peri- 
odic orbit. The usual WKB approximation is the special 
case with two conjugate points v — 2. 
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FIG. 3. The periodic orbit in tlie presence of botlr soft and 
hard confinement potentials. The filled circle is the reflection 
point, and the filled square is the conjugate point. The differ- 
ence between them can be seen quite clearly by plotting the 
coordinate x{t) as a function of time t. 

In the presence of the hard-wall boundaries, we need 
to construct a new quantization rule that includes the ex- 
tra Berry phase. Followi ng a similar calculation, the new 
quantization rule in Eq.(l.l) is derived with the modified 
term rir. Since we have treated the hard-wall boundary 
properly, the modified EBK quantization rule can be ap- 
plied to physical systems with both smooth (soft) and 
sharp (hard) confinement potentials. To illustrate this 
point, let us consider a particle in a harmonic trap, but 
with a hard wall boundary caused by the sample edge as 
shown in Fig. 3, 



V{x) 



2 kx J X ^ 
oo, a: < 



(5.2) 



The classical periodic orbit, shown in Fig. 3, has one 
conjugate point v = 1 and one refiection point r = 1 
which differs from the assumption of the WKB approx- 
imation that there are always two conjugates points. In 
the semiclassical limit, the quantized energy is given by 
the new quantization rule. 



2 ^ 2ni{E - ^kx^)dx = 2tt{ 



+ (5.3) 



where Xt = \/2E/k is the conjugate (turning) point. The 
integral is elementary and leads to the correct energy 
levels 



En = {2n+-)u;, 



(5.4) 



where n = 0, 1, 2, . . . and lo — ^Jkjm is the natural fre- 
quency. The result can of course be verified by the parity 
argument. We emphasize that neither Bohr-Sommerfeld 
quantization (i^ = 0, r = 0) nor WKB approximation 
(i^ = 2, r = 0) leads to the desired result. 

In many practical systems, e.g., quantum dots, quan- 
tum wells, electronic wave guides, etc., the presence of 
both soft and hard potentials is inevitable. The mod- 
ified EBK quantization rule provides us with a conve- 
nient tool to estimate the energy level without solving 
the Schrodinger equation directly. For example, in the 
quantum well biased by a voltage V (that needs not to 
be small) across the well, if the energy is below the bias 



voltage V ^ the periodic orbit has one reflection point and 
one conjugate point. If the energy is above the bias volt- 
age, we have two refiection points. One should not be 
puzzled that the infinitely sharp hard-wall does not ex- 
ist in any practical systems and the correction due to 
hard- wall boundary is only an artifact. As long as the 
potential profile is sharp in comparison to the relevant en- 
ergy gradient of the interested system, the Berry phase tt 
correction is reasonably good. On the other hand, if the 
potential profile is smooth in the same sense, the Berry 
phase should be 7r/2 as in the usual WKB approximation. 

The modified EBK quantization rule can also be ap- 
plied to physical systems in higher dimensions. Let 
us consider a spherical or hemispherical qunatum dot. 
We can either apply the modified EBK formula directly 
to the true three-dimensional trajectories JTst or apply 
the formula after reducing the system to one dimension. 
Here we adapt the second approach. After separation of 
variables, the radial effective Hamiltonian of the three- 
dimensional spherical (hemispherical) quantum dot be- 
comes one-dimensional with the effective potential 



V = 



iji+i) 

OO, r > 



r < 



(5.5) 



where I is the quantized angular momentum. For the 
spherical quantum dot, / takes on all integer values, while 
for the hemispherical dot, only odd integers are allowed 
due to the fiat boundary. 

The classical trajectory of the electron is confined be- 
tween the hard- wall boundary at the surface and the cen- 
trifugal potential near the origin. Thus, there are one 
refiection point r — 1 and one conjugate point ~ I. 
Applying the modified EBK quantization rule, the ap- 
proximate energy satisfies the algebraic equation. 



rE) 



1 



y^O+T)' 



(5.6) 



where rs = + l)/{2mE) is the conjugate point 

and a is the radius of the dot. Instead of solving the 
Schrodinger equation directly, the energy levels can be 
determined easily by the algebraic equation. In the semi- 
classical limit, the conjugate point is close to the origin, 
i.e., a/rE 3> 1. The approximate expression can be fur- 
ther simplified. 



.,.3 1' 
2ma^ V 4 2 



(5.7) 



where V = ^Jl{l + 1). 

Notice that this problem can be solved exactly by the 
spherical Bessel functions. The hard- wall boundary re- 
quires the wave function vanishes at the surface of the 
sphere, ji{V2mEa) = 0, that leads to quantized energy 
levels. In the same limit a/rs ^ 1, the spherical Bessel 
function is approximated by the asymptotic expansion 
that leads to 
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Tpex 



1 2 



n+- 



(5.8) 



The abo ve exact result does not seem to agree with 
Eq. (5.7) at first glance. However, if the angular mo- 
mentum is also semiclassical {I 3> 1), the last term in 
Eq. (113) is ~ 1/2 + 1/4 up to 0{l/l) corrections. It 
is then clear that both give the same result. We empha- 
size again that the agreement is only possible when the 
modified term due to the hard-wall boundary is included. 



Finally, we arrive at the useful Poisson summation for- 
mula, 



/(-«) = -E^(— )■ 



(A3) 



VI. CONCLUSIONS 

In this paper, we study the Berry phase of classical 
trajectories due to the hard wall boundaries. It is shown 
that, upon each reflection by the hard wall, a Berry phase 
TT accumulates. We also relate the hard wall boundary ap- 
proach to the quantum statistics and the discrete gauge 
symmetry. A new quantization rule is derived from the 
modified Van Vleck's formula and applied to simple ex- 
amples. Unlike the WKB approximation that is only 
applicable to smooth potential profiles, the new quanti- 
zation rule provides us with an easy way to estimate the 
energy levels in the presence of both smooth and sharp 
confinement potentials. 
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2112-M-007-101 (HHL), 90-2112-M-007-027 (HHL), and 
89-2112-M007-076 (TMH). 



APPENDIX A: POISSON SUMMATION 
FORMULA 

Poisson summation formula provides a convenient way 
to related two infinite summations together. Let us con- 
sider a physical system on a finite ring with length L and 
lattice constant a. The total number of sites is iV = L/a. 
The discrete version of the usual delta function is 
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^ G=2n7r/a 



(Al) 



where G is the reciprocal lattice vector, 
following summation. 



Akx 



Consider the 



(A2) 



where Xn — na and F[k) is the Fourier transformation 
of f{x). With the help of the identity in Eq.(Al), the 
summation over coordinates is turned into another sum- 
mation over reciprocal momenta. Taking the thermody- 
namical limit L oo, the discrete delta functions are 
related to the continuous ones by LSk,G = 2Trd{k — G). 
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